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AN ANALYSIS OF THE PRACTICAL DPG METHOD 

J. GOPALAKRISHNAN AND W. QIU 

Abstract. We give a complete error analysis of the Discontinuous Petrov Galcrkin (DPG) 
method, accounting for all the approximations made in its practical implementation. Specifi- 
cally, we consider the DPG method that uses a trial space consisting of polynomials of degree 
p on each mesh element. Earlier works showed that there is a "trial-to-test" operator T, 
OQ ' which when applied to the trial space, defines a test space that guarantees stability. In DPG 

J>~y formulations, this operator T is local: it can be applied element- by-element. However, an 

infinite dimensional problem on each mesh clement needed to be solved to apply T. In prac- 
tical computations, T is approximated using polynomials of some degree r > p on each mesh 
element. We show that this approximation maintains optimal convergence rates, provided 
that r > p + N, where N is the space dimension (two or more), for the Laplace equation. 
We also prove a similar result for the DPG method for linear elasticity. Remarks on the 
conditioning of the stiffness matrix in DPG methods arc also included. 
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1. Introduction 



In this paper we prove error estimates for the discontinuous Petrov- Galerkin (DPG) 
method applied to the Laplace equation and the equations of linear elasticity. The approach 

CNJ ■ is applicable more generally to other equations as well. An error analysis of an "ideal" DPG 

method was provided in [5]. Although the ideal method is not practically implementable, 
a number of important theoretical tools for analysis were developed in [5]. We extend this 

CN " analysis using a few new lemmas to provide a complete analysis of the fully implementable 



"practical" DPG method. The distinction between the ideal and practical methods will be 
clear in the next few paragraphs. 

Both methods are easy to describe in a general context. Suppose we want to approximate 
"?/ G U satisfying 

b(v,r) = l(r), WeV. (1.1) 

^ ■ Here U is a Hilbert space with norm || ■ \\u and V is a Hilbert space under an inner product 

(•, -)v with corresponding norm || • \\y. (All spaces are over R.) We assume that the bilinear 
form b(-, •) : U x V i— > R is continuous and the linear form /(■) : V \-t R is also continuous. 
Define T : U H- V by 

(TV, t) v = b(w, y), W E V. (1.2) 

Then, the DPG approximation to w, lies in a finite dimensional trial subspace Uh C U (where 
h denotes a parameter determining the finite dimension). It satisfies 

b(vh,r) = l(r), WeV h , (1.3) 
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2 J. GOPALAKRISHNAN AND W. QIU 

where Vh = T(Uh). Since Uh 7^ Vh i n general, this is a Petrov-Galerkin approximation. The 
method (jl.3p is the ideal DPG method. The excellent stability and approximation properties 
of this method are well known pEJ [5] . 

The main difficulty of the ideal method is that in order to compute "?/h, one needs a basis 
for Vh, which must be obtained by applying T. This is infeasible, as seen from (II. 2p . if 
V is infinite dimensional, unless a solution to ( II. 2p can be written out in closed form. In 
certain one-dimensional problems, and in some multi-dimensional problems like the transport 
equation, the application of T can be exactly written out in closed form (see OH]). But for 
the vast majority of interesting problems, this is not possible. 

Yet, one may approximate T by T r , defined as follows. Let V r be a finite dimensional 
subspace of V (where r is a parameter determining the finite dimension). Then T r w in V r 
is defined by 

(TV, y) v = b{w, y), Vr e V r . (1.4) 

One can then reconsider the DPG method (II. 3p with V h r = T r (Uh) in place of Vh, i.e., 

6KV) = /(r), WeV£. (1.5) 

This yields an implementable method that is very generally applicable. We refer to this 
method as the practical DPG method. 

A serious difficulty still remains when these ideas are applied to standard variational 
problems. Namely, one application of T r requires inverting a Gram matrix in the ^-inner 
product. This is prohibitively expensive for most standard variational formulations. For 
instance, if V is if 1 (i7), where Q is the computational domain, then inverting the Gram 
matrix is as expensive as solving the Laplace's equation. 

This difficulty can be overcome by hybridization, as shown in the earlier DPG papers [HIS]. 
Namely, given a boundary value problem, introducing certain interelement fluxes and traces 
as new unknowns, we can design an ultraweak well-posed variational formulation involving 
a space V that contains functions discontinuous across mesh element interfaces. This then 
implies that the Gram matrix becomes block diagonal, with one block per mesh element 
(since V r may now be chosen to be a DG subspace). The application of T r is thus reduced 
to an easy block diagonal inversion, i.e., the action of the operator T r is local. 

Such an ultraweak variational formulation has been developed for the Poisson equation 
in [5], where its wellposedness is also proved. We will heavily rely on such wellposedness 
results in this paper. An ultraweak formulation for the linear elasticity system is also avail- 
able now [SJ. Both these works analyzed the ideal DPG method (11.31) for the respective 
ultraweak formulations. The aim of the present paper is to provide an error analysis for the 
corresponding practical DPG methods (II. 5p . 

In the next section we will present an approach to the analysis of the practical method, 
continuing in the general context and using the abstract notations introduced above. In 
Section [3], we perform the error analysis for the practical DPG method for the Laplace 
equation. We also provide a condition number estimate. In Section HI we consider a second 
example of linear elasticity and provide an error analysis. 

2. The approach to analysis 

The purpose of this section is to point out a simple functional analytic route to proving the 
discrete stability of the practical DPG method (II. 5p . This discrete stability will follow from a 
discrete inf-sup condition on the space V[. However, in applications, it is often inconvenient 



AN ANALYSIS OF THE PRACTICAL DPG METHOD 3 

to work directly with this space. We prefer to work with V r , which will be some standard 
polynomial space in most applications (as seen in the examples later). The next theorem 
shows that the existence of a Fortin operator into this more standard space V r is a sufficient 
condition for the convergence of the practical DPG method. 

Before we give the result, let us state the assumptions that we shall verify for each of our 
examples. We assume that 

{w 6 U : &(*-, f) = 0, Vr eV} = {0} (2.1) 

and that there is a positive constant C\ such that 

CilMlv < sup h -^P~, Vr G V. (2.2) 

Above and throughout, we will tacitly assume that the suprema such as the above are taken 
over nonzero functions. Let C<2 > be such that 

b(r,y)<C 2 \\w\\u\\-r\\ v , \fweU,reV. (2.3) 

Clearly, such a Ci exists due to the continuity of b(-, •). Finally, assume that there exists a 
linear operator 77 : V H- V r such that for all v G V, we have 

b(w, r - lit) =0, v> G U h , (2.4a) 

\\nr\\ v < C n \\y\\ v . (2.4b) 

Theorem 2.1. Suppose the assumptions (12. ip . (12.21) . (12.31) . and (j2.4p hold. Then the prob- 
lem (11.11) is well-posed and 

\\*-*h\\u< ^ mf \\*-*\\u. (2.5) 

Proof. We apply Babuska's theory [U [TTJ. Accordingly, if we prove the discrete inf-sup 
condition 

^\Mu < sup b -^p-, Ww G U h , (2.6) 

C77 reV£ ||r||y 

then (I2.5P will follow. We prove (I2.6P in three steps, the first two of which are fairly standard 
(but included for readability). 

As the first step, we prove that the following inf-sup condition holds: 

CiM\ v < sup l ^ y ;. y)l , Vr G U. (2.7) 

rev IM|y 

This follows from the other inf-sup condition (12. 2p . Define a linear operator B : U — >• V* 
by Bw = &(#-, •) G V*, for all jr G U. It is well known (12] that (12. 2 p holds if and only if 
B* is injective and the range of B* is closed in U*. Additionally, by (12. ip . B is injective. 
Therefore, by the Closed Range Theorem, B* (V) = U*, so (B*)~~ l : U* — > V exists. Hence 
£-1 . y* _^ jj a j go ex j s ^ s anc [ j s continuous. This proves that problem (11.11) is well-posed. 
We obviously also have WB^W = HCB -1 )*!! = || (J3*)~ x || , i.e., 

-f -p irPP!r = WB-T 1 = IKBTT 1 = mf sup l^j , 

which proves (12. 7\i . 
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As the second step, we prove the following inf-sup condition. 

2-IMI.7 < «P ^1, W <= %. (2.8) 

Note that this differs from (12.61) only in the space in which the supremum is sought. To 
prove (j2.8p . we use (12.71) and assumption (12.41) as follows: 

_ i, b(w,r) biw.FLf) 

CilMli/ < sup -jj— jj — < sup — — 



-rev Mv -rev C n I|i7r| 



Now, since Uv is in V r , the last supremum may be bounded by the supremum over all V r , 
so we obtain (12. SJ) . 

As the third and final step, we prove that if s\ is the supremum in (J2.8P and S2 is the 
supremum in (I2.6p . then si = s 2 - Obviously, si > s 2 as V" 3 V h r . To prove the reverse 
inequality, observe that s\ = \\T r w\\v, by (ll.4p . Since T r w is in V£, we have 

(T r w,T r w) v ^ (T r r,r) v 

si = — — — < sup — r ^ = s 2 . 

||TV||y ?&% \\r\\v 

Therefore, the inf-sup condition (12.61) follows from (I2.8p . □ 

Remark 2.2 (Test basis). The above proof also shows that under the assumptions of Theo- 
rem [2J], the operator T r : Uh M- V r is injective: indeed, if T r w = 0, then b{w, t) = for all 
f in V r , so by the inf-sup condition (I2.8p . we conclude that w — 0. Note that the injectivity 
of T r implies that 

dim(^) = dim(f/ /t ). 
It also implies that a basis for V h r can be computed by applying T T on any basis for Uh- 

Remark 2.3 (Conditioning). Suppose 38 i is a basis for L^. Then, under the assumptions 
of Theorem 12.11 T r ,^ is a basis for V h r , as seen in Remark 12.21 The ij-ih entry of the 
stiffness matrix of the DPG method with respect to this basis is given by Sij = b(,%j, T r ss^) = 
(T r sgj,T r sSi)v- Clearly, S is symmetric. The above mentioned injectivity of T r implies that 
S is also positive definite. To understand the conditioning of S, let us first note that 

T^IMIu < \\ Tr ^\\v < C 2 \\w\\u, Vr e U h . (2.9) 

This follows from the inf-sup condition (I2.8P in the proof of Theorem 12.11 and the continuity 
property f 12 . 3 j) . Next, suppose % = Y^i x i^i is the basis expansion of any % in Uh, and A , Ai 
are positive numbers such that 

Aolkll^ < \\&\\u < -MMI^, Vs" e [4. (2.10) 

Since x t Sx = \\T r x\\u, these estimates imply that the Rayleigh quotient x t Sx/x t x is at most 
AiCf and at least C\\qJC\. Hence 

«(5) < f°g±, (2.11) 

where k(S) is the spectral condition number of S. This gives condition numbers comparable 
to other methods, as we shall see later in our examples. 
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3. First example: Laplace equation 

The ideal DPG method for the Laplace equation was developed and analyzed in [5]. In 
this section, we will set the abstract forms and spaces of the previous section to those from [S] 
and verify the hypotheses required to apply Theorem 12. II Roughly speaking, our main result 
shows that if polynomials of degree p are used to approximate the solution of the Laplace 
equation, then a sufficient condition for optimal convergence is that T is approximated by 
polynomials of degree p + N, where N > 2 is the space dimension. In the wording of j5], this 
means the "enrichment degree" should be chosen to be N. Ample numerical evidence, in 
support of the choice of 2 as enrichment degree, was presented in [HI § 6.1], but all numerical 
experiments were in the two-dimensional case. 

In the remainder of this paper, we let J? be a Lipschitz polyhedron in M. N . We denote by 
{^h}hei a family of conforming shape regular simplicial finite element triangulations of Q. 
The index h now stands for the maximal diameter of simplexes in Q^. 

3.1. Infinite dimensional spaces. Let V = M. N . We use L 2 (^?,V) to denote the set of 
vector-valued functions whose components are square integrable. We set the trial and test 
spaces by 

U = L 2 {Q; V) x L 2 (Q) x Hl ,2 (dQ h ) x H~ l/2 (dQ h ), 
V = H(div,f2 h )xH 1 (Q h ), 

where the "broken" Sobolev spaces (admitting interelement discontinuities) are defined by 
H{div,n h ) = {r : t\ k e H(div,K), MK G fl h } and H\Q h ) = {v : v\ K G H\K), VK G 
Oh}. They have the natural norms 

\\ v \\m(a h ) = (v,v) Qh + (grad^gradw)^, 

IMII(div,/4) = favhh + (divg,divg)^ h . 

The derivatives above, and in such notations throughout, are calculated element by element 
and 

(r,s)a h = Yl (^ s )k, (w,£) 9 a h = ^ (w,£) 1/2 , d K- 
K&n h Ken h 

and (-,£)i/2,dK denotes the action of a functional £ in H~ l l 2 (dK). We will also use \\r\\^ h 
to denote the norm (r,r)^ . The spaces of traces and fluxes are defined by H (dQh) — 
{V e U K Hl/2 ( dK ) :3fflG Hl{Q) such that r)\ 9K = w\ dK VK G L? h }, and H^ 2 (df2 h ) = 
iv e HK^^^idK) : 3q G H(div, ft) such that t)\q K = q ■ n\ dK VK G i?^}, with respective 
norms 



u „i/a 



Hq' (dO h ) 



inf {||iw||jfi(i7) : \/w G Hq((2) such that u\dK = w\qk], (3.1) 

||^n|U-V2(sfi h ) = inf (|M|H(div,fi) : Vg G H(div, Q) such that a n \ 9K = q-n\ dK }- (3.2) 
The spaces U and V are endowed with product norms, i.e., 

\\{a,u,u,a n )\\u= \W\\n + h\\n+ \\u\\ 2 H i/2 {d a h ) + \\vn\\ 2 H -i/2 (dnh) , 
\\(r,v)\\ 2 v = ||T||^ (diViflh) + \\v\\ 2 H i {nh) . 
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3.2. Forms. The ultraweak formulation of the Laplace equation derived in [5] reads as 
follows: Find % = (a, u, u, a n ) G U satisfying (II. ip for y = (r, v) G V where the forms b(-, •) 
and l(-) are set by 

b(it,r) = (a,r)n - (u,divr)a h + {u,r ■ n) dQh - (cr,gradv)n h + (v,a n ) dfJh , 

Kr) = {f,v)a, 

for some / in L 2 (i7). The -u-component of a u solves the Laplace equation with zero Dirichlet 
boundary conditions on df2. For details, consult [5]. 

3.3. Discrete spaces. Let us first establish notation for a few polynomial spaces that we 
will use here and throughout. Let P P (K) denote the space of polynomials of degree at most 
pona simplex K. We write P p (K;Y) for vector valued functions whose components are in 
P p (K). Let A m (K) denote the set of all m-dimensional sub-simplices of K. Define 

P p (K) = {p p G P P (K) : p p \ dK = 0}, 
P p {dK) = {fjL : fi\ F G P p (F), VF G Ajf.xiK)}, 
P p (dK) = P p (dK) n e°(dK), 

where G°(D) denotes the set of continuous functions on any domain D. 

Using these notations, we set the trial approximation space for the DPG method by 

U h = {(a,u,u,a n ) G U: a\ K G P p (K;Y), u\ K G P p (K), 

u\dK G Pp+i(dK), a n \ dK G P p (dK), \/K G Q h }. 

The discrete test space is defined by V£ = T r (Uh), so to complete the prescription of the 
practical DPG method, we only need to specify V r . Set 

V r = {(r, v)EV: t\ k G P r (K; V), v\ K G P r (K), VK G Q h }. (3.3) 

where the degree r > p + N. Clearly, the application of T r , as defined by (ll.4p . can proceed 
locally, element by element, since V r has no interelement continuity constraints. 

3.4. Verification of the assumptions. To apply Theorem I2.ll to the above setting, we 
need to verify its assumptions. 

• Assumption (I2.ip is verified by Lemma 4.1]. 

• Assumption (12. 2p is verified by Theorem 4.2]. 

• Assumption (12.31) is easy to verify. For example, to show the continuity of the term 
(it, t ■ n)da k , we let w G H 1 ^) be any extension of u and observe that 

(u,r-n)dn h = (gradw,r)n h + (w,divr)a h < \\w\\m(n)\\r\\H(div,n h )- 
Taking the infimum over all such extensions w, we obtain 

(u,T-n) d n h < \\u\\ H i/2 {dah) \\T\\ H (div,n h )- 

The other terms in the bilinear form are similar or simpler. 

• Assumption ( 12. 4ft is verified below. 

An operator 77 satisfying ( 12.41) will be constructed in the form Pit = (JJp^r, i7§ rad f). 
We construct the operator i7|p ad in Lemma 13.21 below, and we construct the operator n^ 2 
in Lemma 13.31 But first, we need the following intermediate result. Let Bf raA (K) = {p r G 
P r (K) : p t \e = 0,\/E G Aat_ 2 (A')} and hx = diaxn(K). Hereon we use c and C to denote 
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generic constants (whose value at different occurrences may differ) independent of hx, but 
possibly dependent on the shape regularity of K and the polynomial degree p. We also let 
(•, -)qk denote the L 2 (dK)-miaer product. 

Lemma 3.1. Let r = p + N. Then, for every v G H l (K), there is a unique LT^v G Bf vaA (K) 
satisfying 

{n° r v - v, q p ^) K = 0, Vg p _x G P p ^(K), (3.4a) 

(n?v - v, Hp)dK = 0, V^ p G P p (dK), (3.4b) 

\\H?v\\k + h K \\ gradII?v\\ K < C (\\v\\ K + h K \\ gradv||^) . (3.4c) 

Proof. First, to see that the number of the equations in (I3.4a,p - (13.4b|) equal dimBf rad (K), 
observe that 

dimPf ad (A) = dimP r W + J2 dimp r(F)- (3-5) 

Let bx and bp denote the product of all barycentric coordinates that do not vanish everywhere 
on K and F, resp. Then P r (K) = bxPr-N-iiK) and P r {F) = bFP r -N(F). Therefore, by 
our choice of r, we have dim P r (K) = dimP p _i(A') and dim P r {F) = dim P P (P). It then 
follows from (13. 5p that (I3.4ap - (13.4bj) is a square system for LT®v. 

Hence, to prove that (l3.4ap -f l3T4b|) has a unique solution, it suffices to prove that if v — 0, 
then LT®v = 0. Since LT®v G P r grad (A'), on any face F G A N ^i(K), we may write (LT°v)\ F = 
bpWp for some w p G P P (F). But then, (13.4bj) implies that LT®v must vanish on dK, so 
LT®v = b K z p _i for some z p _ 1 G P p _i(A). Then (I3.4ap implies that LJ^v = on A. 

Finally, one can prove f)3.4cp using a standard affine mapping argument. □ 

Lemma 3.2. Let r = p + N. Define P grad t> = TI®{v — v) +v, where v\k = |A| _1 J K V - Then 

(Lir d v - v, q p ^) K = 0, Vg p _x G P p -i(K), (3.6a) 

(nr d v - v, fi p ) dK = 0, V^ p G P p (dK), (3.6b) 

\\n? ad v\\HHK) < C\\v\\ HHK) , Wv G H\K). (3.6c) 



Proof. Obviously, TIf raA v-v = (TI?-I)(v-v). Hence, flSTBaj) and ( l3~6bl) follows from (HOajl 



and ( 13.4b)) of Lemma 13.11 It remains to prove (I3.6cp . By (I3.4cl) and the Poincare-Friedrichs 
inequality, 

\\nr d v\\K<\\v\\K + \\n°(v-v)\\ K 

< \\v\\ K + C (\\v - v\\ K + h K \\ grad(f - v)\\k) 

< C(\\v\\ K + h K \\gradv\\ K ) , and 
h K \\ grzdnr d v\\ K = h K \\ grad i7> - v)\\ K 

< C r (||v-iJ||j{- + /ijic||grad(i;- v)||a-) 

< Chx\\ grad t; 1 1 a-. 

Canceling out Kk and adding, (I3.6cp follows. D 
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Lemma 3.3. There is an operator P d + 2 : H(div,K) \— > P p+2 (K;Y) such that for every 
t G H(div,K), we have 

(77 d >, q p ) K = (r, q p ) K , \/q p G P p (K; V), (3.7a) 

(P d +V ■ n, Hp+ifQK = (fJ-p+i, r ■ n) 1/2 ,dK V// p+ i G P p+1 (dK), (3.7b) 

H-^p+2' r l|H(div,X) < C||T||H(div,X)- (3.7c) 

Proof. We will first construct the operator on the unit simplex K in R . Recalling the 
notations in §3.31 define P p (dK) to be the L 2 (<9i^)-orthogonal complement of P p (dK) in 
P p (dK), and 

B%{K) = {t e P P +i(#; V) + xP p +i(A) : (p±, f • n>^ = 0, \/p ± G P^*)}. 
We construct an operator P d + 2 mapping H(div, K) into Bp" 2 (K) by 

(r^_ v 2 f,<? P ) A = (r,q P )K, Vg p G P P (K; V), (3.8a) 

{n p + 2 f ■ n, Ap+iW = (Ap+i, r ■ n)i/2,dK V/lp+i G P p+ i(<9#). (3.8b) 

We claim that (I3.8ap - fl3.8bp uniquely determine IJ p " 2 r G -Ey^-^O- Indeed, if their right 
hand sides vanish, then since n p " 2 f is in B p " 2 (K), we find that P p + 2 is a function in the 
Raviart-Thomas space whose canonical degrees of freedom vanish (see e.g., [101 Definition 5]), 
so fl^ 2 r = 0. Hence fl378aD - fl3~8bj) uniquely defines flf^ 2 f. 

Now, we define P p + 2 on any general simplex K by mapping TI p " 2 f from K using the 
Piola transform, as follows. Let Gk be the affine homeomorphism from K onto K and let 
A denote its Frechet derivative. Given any r G H(div, K), let f(x) in ET(div, i^) be defined 
by r o G K = (det A) -1 Ar. Then, define i7 d ; v 2 r by 

PS-Xx) = ^—^n% 2 t{x), with x = G*(s). 

We will now show that this P p + 2 r satisfies the three properties in (13. 7p . 

First, observe that f)3.8al) and f!3.8bp imply the corresponding identities on K, namely, 

(77 d > - r, A-% o G K }) K = 0, Vg p G P P (K; V), 

(n p " 2 r ■ n, fi p+1 o Gf<) dK = (/l p+ i o G K X , t ■ n) 1/2>dK V/x p+ i G Pp+it&K - ). 

This implies ( 13.7a[) and ( I3.7bh . 

It only remains to prove f)3.7cl) . We do this in two steps. First, we prove an L 2 (K) bound 
using the Piola map's well-known estimates for shape regular meshes, namely 

IK] 1 / 2 „ „ 

c\\r\\ k < \\t\\ k — <C\\r\\ k , 

n K 

4 divf||^ < || divr||^|K| 1/2 < C\\ divf ||#. 
Together with the fact that P p + 2 * s a continuous operator on H (div, K), we obtain 

l|P p d >lk + Ml divP p d ; v 2 r|U < C(\\r\\ K + h K \\ divr|U) . (3.9) 

In particular, this proves the L 2 (i^)-bound ||P p + 2 r ll-K' — C|| r llH(div,x)- 
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7 div , 

l P+2> 



Next, we prove a better bound on the divergence norm || div LT^YotWk by showing that 



div(77 p d » = i7 p+1 divr 

where LT p+ i is the L 2 (i^)-orthogonal projection onto P p+ i(K). Indeed, for any u p+ i G 
P p+ i(K), we have, due to (I3.7ap and (I3.7bl) . that 

(div(n^ 2 T),u p+1 ) K = -(n^ 2 r, gradu p+l ) K + {(II*% 2 t) ■ n,u p+1 ) dK 
= -(T,gTadu p+1 ) K + (u p+1 ,T ■ n)x/2,dK 
= (div t,u p+1 ) k ■ 

Hence, 

|| div(i7 p d ; v 2 r)|| A ^ = \\n p+1 div t\\ k < || divr||jf. (3.10) 

Estimates ^M) and (ETTUI) prove fl3T7cD . D 

Now we are ready to apply Theorem 12.11 to obtain a convergence result for the practical 
DPG method for the Laplace equation. 

Theorem 3.4. Let r > p + N. Then the exact and discrete solutions for the DPG method 
for the Laplace's equation, namely % = (o~,u,u,a n ) and °//h = (crh,Uh,Uh,cr nt h), satisfy 



\o~ - o- h \\ L 



2 (fi) + ||' u — Uh\\L 2 (Q) + \\u — ^hW^/2^^^ + \\Cr n — <Tn, fe|| #-1/2(914) 



<C inf [\W - p h \\ L 2(Q) + \\u-w h \\ L 2 {n) 

(.Ph,w h ,w h ,fi h )<£U h 



U ~ W h \\ H i/2 + \\a n - Vn,h\\H-i/2(dQ h ) 



1 d 



Proof. As already observed, we have verified the first three assumptions of Theorem 12.11 To 
verify Assumption (12.41) . let r = (t,v) and set LTf = (TIp" 2 T,n^ rad v). The continuity esti- 
mates of n^ 2 and TIf™ A of Lemmas O and E31 (namely fl3^cj) and (ET7cl) ) show that (LTibl 
holds. To see that (I2.4ap also holds, observe that the identities of these lemmas also imply 

(p h , r - n%r) n = 0, K, div(r - n^ 2 r)) Qh = 0, 

(w h , (r - Tlf^r) ■ n) dnh = 0, (p h , grad( W - 77f ad t;)K = 0, 

{v-nr d v,v h ) dnh = o, 

for all (p h , Wh, Wh, f)h) G Uh- While the identities above on the left follow from the identities 
of (13. 6p and (13.71) . those on the right are proved by integration by parts. Together these 
identities imply that b{w, r — Fly) = for all w G Uh, so Assumption (12.41) is satisfied. □ 

Remark 3.5 (Enrichment degree). The above arguments point to the potential of choosing 
different enrichment degrees for the scalar and flux components of the test space. We have 
in fact proved that if, in place of the V r set in (13.31) . we revise our choice of V r to 

V r = {(r,v) eV: r\ K e P P+2 (K;¥), v\ K e P P+N (K), VK g Q h }, 

then, we obtain the same convergence result. Obviously, the revised V r defines a smaller 
space if A r > 3. The present DPG software packages are set to approximate all components 
of T by polynomials of the same degree r. Our results indicate that this is unnecessary. 

As an example of how Theorem 13.41 implies /i-convergence rates, we state the following. 
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Corollary 3.6 (Convergence rates). Let h = max^-g^ diam.(K), N = 2 or 3, and let the 
assumptions of Theorem 3.J\ hold. Then 



11°" — o~h\\L 2 {Q) + \\ u — u h\\h 2 {n) + \\ u — Uh\\ H V2, dn . + \\a n — o- n ^ h \\ H -i/2( dQh - ) 

< Ch s (\\u\\ H s+i( Q ) + HcrHijs+i^)), 
for all 1/2 <s<p + l. 

Proof. The proof proceeds by bounding the infimum over (ph,Wh,Wh,i]h) £ U^ in Theo- 
rem 13.41 It is standard to bound the first two terms in the infimum, so we will only explain 
how to bound the next two terms. It is well-known (see, e.g., [8j Theorem 8.1]) that there 
are interpolants ilgrad^ £ Hq(Q) and Tl^a G H(div,fl), such that TI gTa< iu\K £ P P +i(K), 
ridivO'lK G xPp(K) + P P (K) 3 for all K G O^, and the interpolation errors satisfy 

\\u- n gxSLd u\\ H i^Q) < Ch s \u\ H s+i {n) , (1/2 < s <p + l), (3.11a) 

||o--i7 div o-||H(div,fi) < Ch s \cr\ H s+i(n), (0 < s <p+l). (3.11b) 

Let tft denote the trace of n graA u on dQh- Then, 

inf ||u - w h \\ H i/2 idfih) < \\u - i h \\ H i/2 {dQh) < \\u - n grad u\\ H i {Q) . 

The last inequality is obtained by observing that u is the trace of u on dfih and bounding the 
infimum in definition (13. ip . In a similar fashion, we can estimate the last term in Theorem 13 .41 
by || a — iIdivC"||_ff(div,i?)- The interpolation error estimates (13.111) then finish the proof. □ 

To conclude this section, we prove that the condition number of the stiffness matrix of the 
DPG method is no worse than other standard methods - see Remark 12.31 for the definition 
of the stiffness matrix with respect to a basis {.^}. Consider, for definiteness, the three- 
dimensional tetrahedral case. We tacitly assume that the basis functions &i are local, and 
obtained, as in usual finite element practice, by mapping from the (reference) unit simplex. 
For example, a basis for the trial space for the numerical traces is built using a local basis 
{ej} for P p+ i(dK), which in turn is obtained by mapping over a basis {ij} for P p+ i(dK) 
(where ej = ej o Gj^ and we use the other mapping notations in the proof of Lemma [3. 3p . 
Consequently, if s = J2j s j^j * s the basis expansion for any § G P p+ i(dK), then by the 
equivalence of norms in finite dimensional spaces 

c$>/< inf ||g||* jj <Cj>/. (3.12) 

Such arguments will be used in the following proof without further explanation. 

Theorem 3.7 (Conditioning). Suppose Q^ is a quasiuniform tetrahedral mesh and the as- 
sumptions of Theorem \3.Ji\ hold. Then the spectral condition number of the stiffness matrix 
S of the DPG method satisfies 

k(S) < Ch~ 2 . 

Proof. Let us apply (12. lip . We have already shown above that C\, C 2 and C n are independent 
of h. Hence it only suffices to find the dependence of Aq and Ai on h in (12.101) . 
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Let x = (p,w,z,fj) in U^. As a first step to bound the norm of z, we recall that the 
existence of an H 1 (K) polynomial extension [7] implies that for any S in P p+ i(dK), 

inf ||eJ| 2 > + II gradejl^ < C inf ||e|| 2 > + || grade|| 2 >. 

e P eP p +i(X), eCH^K), 

Mapping to K using the affine homeomorphism Gk and scaling both sides by \K\, we obtain 
inf \\e p \\ 2 K + h 2 K \\ grade p |||- < C inf HeH^ + h K \\ grade||^. (3.13) 

e p eP p+ i(K), eeH 1 ^), 

e P \dK=s e\dK=s 

where 5 = 6*0 Gj^. Let us denote the function which achieves the left infimum by E!= rad s. 
Applying the above inequality, element by element, with s replaced by z, and using Kk < 
diam i?, we have proved that 

\\E^zf n < L7max(l,diam^) 2 inf (||e|| 2 fi + || grade© < C\\z\\ 2 1/2 

(e-2)|ejf=0 



Thus 



4e?*>h\i < piiu < nEf^r 



1 /2 

where the upper inequality is obvious from the definition of the H (<9J?h)-norm. By an 
inverse inequality, 

c\\Ef^z\\l < \\z\\ 2 Hl n {dQh) < Chr 2 \\Ef^z\\l. (3.14) 

A similar argument using the H(div,K) polynomial extension in [8], gives 

c\\E} w f)\\l < \\f}f H -m {d a h ) < Ch^WE^Wl (3.15) 

Combining (13. 14ft and (13.150 . we have 

c|M| 2 < IkH 2 / < Ch- 2 \\x\\l \fx e U h , (3.16) 

where \\x\\ 2 = \\ P f n + \\wf n + \\Ef*zf a + \\E^fj\\l. 

To prove (I2.10p . consider the coefficients in the basis expansion of x. If z/s denote the 
coefficients in a basis expansion of the z\qk, then using (I3.12p and the minimization property 
of E& z, we obtain 

C E N 2 ^ W\ (ll^^llx + h 2 K \\ grad^iliy < C£ | % | 2 . 



By an inverse inequality 

\K 



C / ,1^1 S T^T H-^p Z||jjT S^; y |-2j| • 



A similar estimate holds for E dlY, f). Combining these with the obvious estimates for the 
coefficients in the expansion of p and w, we find that 

c||a;|| 2 2 min \K\ < ||.r||n < C||a;|| 2 2 max \K\. (3.17) 

Clearly, inequalities (13 .17ft and (13. 161) imply (12. 101) with Ao = cmm.Ken h \K\ and Ai = 

Ch~ 2 msLXK£n h \K\, thus completing the proof. D 
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4. Second example: Linear elasticity 

Two ideal DPG methods for the linear elasticity equation were developed and analyzed 
in [2|. The two methods are equivalent for homogeneous isotropic materials. Among their 
many interesting properties is their robustness with respect to the Poisson ratio, i.e., the 
method is locking-free. In this section, we will consider the practical version of one of these 
two methods and prove its optimal convergence. We proceed as in the previous example, by 
first setting the abstract forms and spaces to those specific to this method then proceed to 
verify the hypotheses required to apply Theorem 12.11 In this section, we restrict to N = 2 
or 3. The results and the analysis are similar to those in Section El so we will be brief. 

4.1. The spaces. We set the trial and test spaces by 

U = L 2 (J?;M) x L 2 {Q- V) x H l Q l2 (dQ h - V) x H~ l/2 (dQ h ;Y) x R, 
V = #(div, Q h ; S) x H x (Q h \ V) x L 2 (Q; K) x R, 

where M = M. NxN , § consists of symmetric matrices in M, and IK consists of skew-symmetric 
matrices in M. The trial and test spaces are normed by 

\\(a,u,u,a n ,a)\\l = \\cr\\ 2 n + \\u\\ 2 n + ||u||#i/2 (sr4) + \\^n\\ 2 H -i/2 {mh) + \a\ 2 , 
\\(T,v,q,P)\\lr = ||r||^ (divA) + ||u|||n (flh) + \\q\\ 2 n + \f3\ 2 . 



4.2. The forms. The (second) ultraweak formulation derived in [2] reads as follows: Find 

>% = (a, u, it,, a n , a) G U satisfying fll.ip for all v = (r, v, q, (3) G V where the forms &(■, ■) and 
/(■) are set by 

b(w,-y) = (Aa,r)n h + (u,divr)n h - (u,rn) dnh +Q 1 (aI,AT) n 
+ 0, gradv)n h + (o", q)n ~ («, v n )dn h 
+ Q 1 (Aa,/3I) Q 

l(*) = (f,v) n , 

for some / in L 2 (i7;V). Here, A is the generalized compliance tensor (see e.g., [21 Re- 
mark 2.1]) and Qq is the essential infimum of the trace of the matrix A(x)I over x G Q. 
Throughout, we assume that A is element- wise constant. We note that above and through- 
out, the inner products of matrix- valued functions, such as (a, r) K , are computed by inte- 
grating the Frobenius product of the two matrices. 

It is easy to see that the resulting a and u satisfies Aa = e(u), where e(u) = (gradn + 
(gradu)')/2, and diver — f on f2, and u = on <9l?, and a = 0. For details, consult [2]. 

4.3. Discrete spaces. Symmetric, skew- symmetric, and general matrix-valued functions 
whose entries are in P p (K) are denoted by P p (K; §), P p (K; K), and P p (K;M), resp. Using 
these notations, we set the trial approximation space for the DPG method by 

U h = {{a,u,u,a n ,a) G U: a\ K G P p (K;M), u\ K G P P (K;Y), 

u\dK G P P+ i(dK; V), o n \ m G P p (dK; V), «6l, \/K G Q h }. 
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The discrete test space is defined by V[ = T r (Uh), so to complete the prescription of the 
practical DPG method, we only need to specify V r . Set 

V r = {(r,v,q,(3) G V : r\ K G P r {K;S), v\ K G P r (K;V), 

q\ K G P P (K; K), /3 6 1, VK G !l h }, 

for some integer r > p + N. 

4.4. Verification of the assumptions. To apply Theorem 12.11 to the above setting, we 
need to verify its assumptions. 

• Assumption ( 12. ip is verified by [21 Lemma 5.2]. 

• Assumption ( 12. 2ft is verified by [21 Lemma 5.3]. 

• Assumption ( 12.31) can be easily verified, as in the case of the Laplace equation. 

• Assumption (j2.4[) is verified next. 

Let f = (r, v, q, /3) G V. The operator 77 satisfying (12.41) will take the form 

/Jr = {n<$*>r, n^v, tifq, (3). (4.1) 

We set n& Tad v to be the one defined in Lemma |3~2| but applied component by component, to 
the vector valued function v. The operator Z7? is simply the L 2 -orthogonal projection onto 

{q G L 2 (f2;K) : q\ K G P p (K; IK), \/K G i?^}. It remains to construct the operator 11^^' . 
We do so, based on a set of degrees of freedom given in [9], in the next lemma. 

Lemma 4.1. There is an operator II ^ 2 '■ H(div,K]S) — > P P +2(K;E>) such that for every 
t G H(div, K; §), we have 

(i7^2' S) r, q p ) K = (r, q p ) K , Vq p G P p (K; S), (4.2a) 

(i7^2' S) r • n, fi p+ i) d K = (fip+i, r ■ n)x/2,8K, tytp+i G P p+ i(dK; V), (4.2b) 

II j-j-(div,S) || ^ ^-(|| || f A n \ 

\\ U p+2 T\\ H ( diVjK ) <C\\T\\ H ( d i VjK ). (4.2c) 

Proof. We only give the proof for N = 3 as the proof for N = 2 is similar. As in the proof 
of Lemma |3.3[ we will first construct the operator on the unit simplex K in M. N . Define 
P^- +1 (dK;Y) = L 2 (dK; V)-orthogonal complement of P p+1 (dK;Y) in P p+1 (dK;Y) and set 
P° +2 (K;S) = {f G P P+2 (K;S) : (s,frT • n+)e = 0, Vs G P P+2 (e), Ve G Ai(£)} where, 
for each edge e G Ai(if), n + and h~ are the normal vectors of the two faces sharing e. 
Let B^ 2 {K;S) = {f G P p ° +2 (Jl";§) : {v,rh) ak = for all u G P^dK; V)}. We define 

77^2' S) = ff(div, £;S) H- B$ 2 (K;S) by 

(i7^ 2 v ' s) f, g p )^ = (f, g,)*, Vg p G P p (£; S), (4.3a) 

^div,§) . _ ^ fi p+1 ) dk = (/ip+i, f ■ n)x/ 2 ,dK, V/tp+i e P p+ i(dK; V). (4.3b) 

By [3 Theorem 2.1], these equations are uniquely solvable, so Pp+2' is well defined. 

Next, we define -^+2' on an y general simplex il" by mapping II p ^' from il" using the 
Piola transform for symmetric matrix-valued functions. Recalling the mapping Gk from K 
onto K and its derivative A, we define 
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for any r G H(div, K\ §). Here, given r on K, the function f on K is defined by (det A)t(x) = 
Af(x)A t , with x = Gk(x). As in the proof of Lemma T3.3[ it is now easy to see that II Jj+J r 
satisfies (l4T2a) and (OB . 

Next, we observe that the commutativity property 

div il^V = n p+1 div r, (4.4) 

holds, where II p+1 denotes the L 2 (K; V)-orthogonal projection onto P p+ i(K, V). Let w p+ i G 

P p+1 (/l;V). Then 

(div(i7^2' S) r),a;p +1 )^ = -(il^'^T, gradw p+1 ) A + ((i^'^r) • n, u p+1 } dK 

= -(^p+2 ,S) ^ £ (^p+i))k + ((-H^S^M ■ n,W p+ i) 9A ' 
= -(r, e(cj p+ i)) x + (wp+i, 7- ■ ri)i/2,dK, (by (EI}), 
= -(r,gradwp + i)i<: + (u p+1 ,t ■ n) 1/2idK , 
= (divr,uj p+1 ) K . 

which proves (I4.4p . 

It only remains to prove the estimate of (I4.2cp . This can now be done as in the proof of 
the estimate (I3.7cp of Lemma 13.31 in two steps, using (I4.4p in place of the commutativity 
property used there. □ 

The main result of this section is the following. 

Theorem 4.2. Suppose that r > p + N and suppose that the compliance tensor A is 
element-wise constant. Then, the difference between the discrete solution of the practical 
DPG method, %^ = (o , h,Uf l ,,Uf l ,,o' nj f l ,ah) ) and the exact solution a u = (a,u,u,a n ,a) satisfies 

||°" _ 0~h\\L 2 {n) + \\u — Uh\\L 2 {Q) + ||W — ' a h\\ H ^/2, dQ , + \\a n — <7n,7i II #-1/2(514) + l a ~ a h\ 
<C inf I || Cr - Ph\\l?{Q) + \\U - W h \\ L 2(Q) 

(ph,w h ,w h ,fi h ,)eu h \ 

+ ||« - «fcHjfV a (8fl h ) + \\0~n - ffn,fc||ff-i/3(8rt fc ) 

Proof. As mentioned above, we only need to verify Assumption (12.41) for the II in (14. ip and 
apply Theorem 12.11 By the inequalities of previous lemmas, the estimate (12.4b|) is obvious. 



To prove (I2.4al) . namely b(w, f — II y) = for all w G Uh, it suffices to prove the following 
eight identities 

(Ap h , t - n^l' s 'r)a = 0, (w h , div(r - n^l' s) r))n h = 0, 

(w h , (r - n^l'%) ■ n) dQh = 0, (p h , grad( W - II^v))^ = 0, 

(v - nr d v, f) h ) d n h = 0, Qo'ilhl, A(r - n^r))o = 

(p h , q - nfq) n = 0, Qo\Ap h , (P ~ P)I)t> = 0, 

for all w = (ph, w h,ibh,Vh,7h) £ Uh- The first five are proved exactly as in the proof of 
Theorem 13.41 but using the new lemma. The sixth is obvious from (I4.2ap . To see the 
seventh, denoting by skwp^ the skew-symmetric part of ph, observe that (ph,q — n^q)^ = 
(skw(p h ), q - nfq) n = 0, by the definition of Ilf. 
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Therefore, applying Theorem 12. 1[ we obtain a quasioptimality estimate. This yields the 
estimate of the theorem after observing that in the infimum over w = (ph, Wh,Wh,fih, jh) i n 
Uh, we may choose 7^ = a. □ 

We conclude by noting that results similar to Corollary 13.61 and Theorem 13.71 can be 
established for this example as well. The arguments are very similar. 
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